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of fluorescence in the bacterial reaction center induced by entangled light demonstrate a degree of selectivity of 
double-exciton states which is not possible using classical stochastic light with the same power spectrum. 
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I. INTRODUCTION 

Apart from their evident importance in experimental tests 
of the foundations of quantum mechanics, entangled photons 
promise many applications to quantum information process- 
ing Q3I21, secure quantum communication BHD, lithography 
(5HH or metrology fT0HT2ll . In addition, their non-classical 
frequency and time correlations could also open up novel 
spectroscopic applications (T3J PHI . Entanglement-induced 
two-photon transparency fl5\ and the linear scaling of two- 
photon induced fluorescence with the pump intensity lfT6l 
constitute two basic non-classical effects observed with entan- 
gled photon pairs. This scaling makes it possible to carry out 
nonlinear optical measurements with much lower light inten- 
sity compared to classical light, as we will discuss in section 
|lllj More generally, entangled photons offer new control pa- 
rameters for nonlinear spectroscopy [17], and can be used to 
distinguish quantum pathways of matter fl8l . It remains how- 
ever an open question, to what extent quantum entanglement 
is essential for these effects and whether some can be repro- 
duced, for instance, by shaped or stochastic classical pulses 
Il9l|20]l. Recently, it has been argued that entanglement can 
induce collective resonances between non-interacting two- 
level atoms I2TII . It was later found that these resonances may 
not be observed in fluorescence due to destructive interference 
of pathways fT3ll . 

In this Paper we present a theoretical framework for describ- 
ing excitations induced by entangled photons in chromophore 
aggregates. These possess manifolds of single- and double- 
exciton states 1 22 ], and the excitation of each double-exciton 
state can take place via many interfering pathways (see figure 
[TJ. We show that the non-classical spectral and temporal prop- 
erties of twin photons can be used to excite state distributions 
that may not be obtained by classical light. Due to energy con- 
servation in their generation, the sum of the energies of twin 
photons u p = uj\ + UJ2 has a sharp distribution (i.e. Au p is 
very small), even though each single photon may be broad- 
band (Acji, Aco>2 very large). The high frequency resolution 
of uj p allows to target certain double-exciton states, while the 
broad bandwidth of the individual beams enables the simulta- 
neous access of all possible pathways that reach the double- 
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FIG. 1. (Color online) A: Two-photon transitions with entangled 
photons. The broad bandwidth allows us to access the entire man- 
ifold |e), at the same time the manifold |/) is well-resolved due to 
the narrow distribution of uj\ +CJ2. B: Transitions with classical light. 
Bandwidths of uj\ and 002 now add up, and the f-manifold cannot be 
resolved. 



exciton state (see figure [T}. In contrast, with classical fields 
the bandwidths of the two beams add up, and create a large 
uncertainty in the sum. Our simulations show that the non- 
classical distributions could be detected by the fluorescence 
signal. The control parameters of the entangled light can be 
used to stretch the signal in two-dimensional plots, that could 
reveal additional information about the matter. 



II. THE MODEL 

We consider a chromophore aggregate interacting with the 
electromagnetic field and described by the Hamiltonian 
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Hq, Hp and H[ nt represent the aggregate, the field and their 
coupling, respectively. The electronic states group into well- 
separated manifolds, which are denoted as e- (single-exciton), 
f- (double-exciton) manifold and so on. The field couples to 
the system via the dipole operator, which induces transitions 
between these manifolds. In the rotating wave approximation, 
the interaction Hamiltonian reads 



H M {t) = V(t)Ei{t) + V\t)E{t), 
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where we have introduced the positive-frequency component 
of the dipole operator 
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in which ujij — (Ei—Ej)/K are matter transition frequencies, 
and fiij the dipole moments. The corresponding negative- 
frequency part of the electromagnetic field operator is given 
by 
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where Q denotes the quantization volume, a\ is the creation 
operator for mode s, and s runs over the relevant modes. As- 
suming perfect phase matching, we can neglect the spatial de- 
pendencies of the electromagnetic field. This is justified for 
large samples. Coupling with a phonon bath can induce many 
interesting relaxation effects. These are not included in the 
present study. We first derive general expressions for density 
matrix elements in the single- and the double-exciton mani- 
fold in terms of convolutions of matter- and field-correlation 
functions. These results are valid for light with arbitrary sta- 
tistical properties, and can describe two-photon induced flu- 
orescence and related measurements in a consistent manner. 
The matrix element (e^, ej) of the single-exciton part of the 
density matrix is given by 



.(t)=tr[\e j (t))(e i (t)\g(t)], 



(5) 



where g(t) denotes the density matrix of the entire matter- 
field system. In a superoperator Liouville space representa- 
tion, the formal solution of the Heisenberg equation is given 
by the Dyson series l23l 
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where T is the time-ordering operator, and the superoperator 
#i nt) _ is defined through A_g = Ag — gA in the interac- 
tion picture with respect to Hq + Hp. For weak fields we 
can expand the exponential, and the leading (second-) order 
contribution to eq. ^ yields 
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where A^ (t) = | ej (t) ) (e^(t) |, and T denotes the set of control 
parameters that govern the excitation, e.g. frequencies, pulse 
envelopes etc. The square brackets (...) denote the quantum 
mechanical expectation value with respect to the initial state 

g(-oo), 



g(-oo) = \g)(g\ <g> ^ fi eid, 
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i.e. the aggregate is in the ground state \g), and the electro- 
magnetic field is prepared in a state £>fi e id. 
The fourth-order processes that end in the \e) manifold can be 
represented by the diagrams I and II shown in figure [2] (for di- 
agram rules, see l24l ). These can be divided into two classes: 
Raman type that only involve single-exciton resonances (dia- 
gram I), and two-photon absorption type (diagram II) that in- 
volve the double-exciton manifold. We have to take two addi- 
tional diagrams into account, in which three interactions occur 
on the right-hand side. Apart from the projector | e j (t) ) (e^ (t) | , 
these are simply the complex conjugates of diagrams I and II 
in figure [2] We can thus write them as 
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FIG. 2. (Color online) Diagrams I and II: the single-exciton density matrix. Diagram III: the two-exciton density matrix. 



e.-, we can combine the two terms 



Note that for diagonal density matrix elements, i.e. 

as the real part of the first term. To the same order, double-exciton matrix elements are given 
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;n by diagram III of figure [2] 
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where B^ (i) 



\fj(t))(fi(t)\. Eqs. ([9), ([TO]) and ([TTJ enable 
us to evaluate the density matrix for arbitrary states of the ra- 
diation field once we specify its four-point field correlation 
function. 



III. LIGHT SOURCES 



A. The twin photon state 



Extensive research effort has focused on producing and de- 
tecting pairs of entangled photons. Numerous schemes have 
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been employed for the generation of entangled photon pairs, 
including parametric down-conversion l25l . biexciton decay 
in semiconductors l26l [27l or four wave mixing in optical 
fibers (281 El. 

We consider a type of entangled light known as twin photons 
(30). It is created by type-II parametric down-conversion with 
a cw-pump laser of frequency u p . When the pump beam is 
sufficiently weak, it predominantly produces temporally non- 
overlapping pairs of entangled photons. The state of the field 
can then be obtained pertubatively in the interaction Hamilto- 
nian of the pump beam with the nonlinear crystal l30l 



\ip) = E p Jdfa Jdfa sinc((/c p - k lz - k 2z )L/2)e i ^- klz - k2z)L/2 

8 (fax + fax) 8 (fay + fay) 8((jj(fa) +U)(fa) - LJ p )\fafa), (12) 



where E p is the pump field amplitude. The state is not nor- 
malized. Here, the length of the crystal along the pump beam 
is denoted L, while it is assumed to be infinitely large in the 
other two directions. In the setup depicted in figure [3] one 
collects the light in two outgoing directions 0\ and 2 with re- 
spect to the optical axis, thus fixing the central frequencies uj\ 



and uj 2 of the photon wavepackets in directions 1 and 2. The 
photon pairs are furthermore correlated in time. This is char- 
acterized by the entanglement time T, which is determined 
by the length L along the axis of the pump beam. Therefore, 
one can vary independently T and two of the three frequen- 
cies co> p , uji and uj 2 . 
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Here, we had defined the rectangular function 



FIG. 3. (Color online) Down-conversion setup to produce twin pho- 
tons: A high-energy pump of frequency u p creates two beams of 
entangled photons with central frequencies cji and CJ2, such that 
Up = cji +CJ2- The central frequencies are determined by the relative 
angles between the direction of propagation. 
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FIG. 4. (Color online) Power spectrum of entangled light [eq. |T3] >] 
with T — 30/s. The central frequencies are cji = 13000 cm -1 
and uji — 11000 cm -1 . The two sine 2 -functions overlap to form a 
single broad peak. 



For short T, eq. ( [T2| ) leads to the broad power spectrum 
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which is depicted in figure [4] To evaluate the normally or- 
dered correlation functions in eq. (10), we note that it can be 
separated into 

(Et(r 3 )S + (r4)£7(r 2 )£7(ri)) = 
(^E^^E^T^iOlE^EiT^), (14) 

where 130] 
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The non-classicality of the quantum state produced by this 
setup has been verified in numerous experiments I3TH381 . 
One important feature that follows from eqs. (14) and ( p?) 
is that the four-point correlation function of the fields scales 
linearly with the pump intensity \E P \ 2 , while classically it 
scales quadratically l39l . This scaling behavior can be uti- 
lized in experiments to determine intensity regimes, in which 
eq. ( [T2| ) can account for the produced signal. In practice, 
we expect this linear scaling of fourth-order signals with a 
crossover to a quadratic classical scaling at higher intensi- 
ties. The light beam can be considered to be made of pairs 
of photons. At weak intensities the pairs are temporally well 
separated, and the process is induced by two photons of the 
same pair. At higher intensities it becomes statistically more 
plausible for the two photons to come from different pairs, 
which are not entangled and the classical scaling is recov- 
ered. This crossover has been demonstrated experimentally 
135H381 . For instance, in lf35l the authors detect it at an in- 
tensity (of the entangled photon flux) of 1.5 fiW in a PPKTP 
crystal, which was obtained by a pump intensity of less than 
2.5 W. The linear scaling of the two-photon absorption was 
also demonstrated experimentally for organic porphyrin den- 
drimers £23133. 

Entangled photon sources are weak, but at the same time weak 
intensities are required to see the effect of entanglement. Non- 
linear spectroscopy and imaging with weak fields should al- 
low to limit the damage to biological samples. 



B. Stochastic light with the same power spectrum 

To highlight the effects of entanglement, we consider a clas- 
sical reference state of the field with the same power spectrum. 
To that end, we introduce the frequency decomposition of a 
classical field 



Ei(t) 



(17) 



where Ai and fa are real functions representing the amplitude 
and the phase of the various modes. To reproduce the power 



spectrum eq. (13 ), we took 



(15) 



Ai(cj) = A sinc[(uj-uj 1 )T/2}e- i ^ UJ \ (18) 
A 2 (oj) = A sinc[(uj - uj 2 )T/2]e- iWuj ^ n/2) . (19) 

With this choice <j>{oS) does not affect the power spectrum 
and will be chosen to obtain stationary stochastic light (rather 
than a pulse). We assume that fauj) is a random function, 
so there is no well-defined phase relation between the differ- 
ent frequencies. The quantum mechanical expectation value 
in correlation functions such as (E^ \r^)E^ ' {t^)E{t2)E{ti)) 
then needs to be replaced by an average over the distribution 



5 



A B C 



special pair 

CO 




FIG. 5. (Color online) A: Arrangement of the six chromophores of the bacterial RC (PDB code: 1PRC) |40|. B: The absorption spectrum of 
the RC. The lifetime broadening of the excited states is taken to be 200 cm -1 . The gray, dashed line shows the power spectrum of figure|4] the 
single exciton states are indicated by red lines, and their heights are proportional to their transition dipole moments \iJL ge \ 2 • C: Level scheme 
of the RC. The 12 single-exciton states are marked in red, the 41 double-exciton states in blue. 



of (j)(uj): 



is one crucial difference. The quantum mechanical expecta- 
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Assuming that <fi is uniformly distributed in the interval [0, 2tt) 
we find that E is a stationary Gaussian process 

(^ t (r 3 )Et(r 4 )^(r 2 )^(r 1 )) 
= <St(r 3 )E(r 2 )><#(r4)^(r 1 )> 
+ (^ t (T4)^(r 2 ))(^ t (r 3 )^(r 1 )). (21) 

Classical light only shows correlations on the level of intensi- 
ties, i.e. (E^E). By substitution of eqs. ( 18 ) and ( 19 ) we then 
obtain 
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(22) 



where the summations run over uj\ and cj 2 , and tri(x) is the tri- 



angular function as defined in (A2). This correlation function 
represents stationary, stochastic light, since it only depends 
on time differences, just like the one produced by entangled 
photons with stationary pump [see eq.(p3)]. However, there 



tion value can be factorized in eq. ( [14] ), and only depends on 
r 3 — T4 and r 2 — Ti . The correlation function ( [22] ) depends on 

ti - t 3 , r 2 - r 4 , n - r 4 and r 2 - r 3 . 



IV. SINGLE AND DOUBLE-EXCITON DENSITY 
MATRICES IN THE BACTERIAL REACTION CENTER 

We now present simulations of the exciton density matrices 
produced by entangled and stochastic light with the same 
spectral density in a model of the bacterial reaction center 
(RC) of purple bacterium B. viridis BH . The primary steps 
in photosynthesis involve excitation energy transfer towards 
and charge separation within RCs. The RC of purple bacteria 
utilizes the high intensity near-IR region of the solar irradi- 
ation, and transforms the energy into a chemical potential 
gradient with near unity efficiency. In bacterial RC the charge 
separation occurs in the active branch of the protein on the 
sub-ps to low ps timescale (42l 03|. The quenching of the 
excited state special pair population has been demonstrated 
experimentally in time resolved emission measurements 
(441143. 



The electronic Hamiltonian Ho describes the optically 
bright chromophore excitations and dark charge separated 
states in the active branch of the RC in a tight-binding 
formulation. It includes 12 single- and 41 double-exciton 
states as depicted in figure [5p. To construct Hq, we start with 
the X-ray structural data for B. viridis (PDB code: 1PRC) 
l40l . The Q y transition dipole moments are placed at the 
center of bacteriochlorophylls and bacteriopheophy tines. 
The excitation energies of BCl's and BP's are taken from 
Ref. 1461 , Forster couplings are calculated in the dipole 
approximation, except for the special pair where a value of 
852 cm -1 was chosen, allowing to reproduce the CD spectra 
of B. viridis. The energy of the primary charge separated (CS) 



FIG. 6. (Color online) Left column: the double-exciton mani- 
fold density matrix upon excitation by entangled light with uj p = 
22160 cm" 1 (A), and with 24200 cm" 1 (B). Right column: same 
but for stochastic light. Only the dominant states /io • • • /30 are 
shown. 



state is fixed relative to the lowest special pair exciton state 
E21 (48), which yields a reference point for the additional CS 
energies of the active branch. The Hamiltonian is given in de- 
tail in ref. (4T].Dephasing effects and incoherent population 
transport were neglected. To model the absorption spectrum 
shown in figure |5j we assume the same lifetime broadening 
7 = 200 cm -1 for all single-exciton states, which agrees 
reasonably with more detailed simulations BTIl . In figure |5p, 
we mark the single-exciton states e by vertical lines, whose 
length is proportional to their oscillator strength |/i^ e | 2 - 
The spectrum is dominated by six states corresponding to 
molecular Frenkel excitations of the constituents of the RC. 
The other six charge transfer states are dark, and cannot be 
accessed spectroscopically. 



FIG. 7. (Color online) Same as figure [6] but for the single-exciton 
manifold. Onlye dominant states . . . eio are shown. 

/30 are most strongly excited. This selectivity allows to ma- 
nipulate of the fluorescence signal, as will be shown in the 
next section. While the density matrices in figures [6] A & B 
strongly depend the pump frequency, our calculations of the 
corresponding density matrices produced by stochastic light 
indicate no such selectivity. The results are plotted in figures 
|6]C & D. For both pump frequencies (ui + U2 = 22160 and 
24200 cm" 1 ), the excitation is distributed among the bright 
states /n, /i 5 , / 22 , /25, f27 and / 30 . 

Figur e [7] shows the single-exciton manifold density matrices 
[eq. ( |A5| )] for the same pump frequencies 22160 and 24200 
cm" 1 . Our results qualitatively resemble the ones for the 
double-exciton manifold. By varying uo p we can enhance (fig- 
ure [7]A) or suppress (figure [7^) the excitation of e$ (and, vice 
versa, of e§ and ej), while stochastic light always creates very 
similar density matrices. For both frequencies, the population 
is mainly distributed between the states e$ and ej. 



The density matrix induced by entangled or stochastic light 
is block-diagonal in the e- and f-manifolds, since a coherence 
between the two manifolds can be related to field correlation 
functions of the kind (E^E^E) or (E^EE), which vanish for 
any Fock state or stationary Gaussian process. In figure [6] we 
depict the absolute values of the density matrix elements of the 
double-exciton states fio to /30 for uj p = 22160 cm" 1 and 
24200 cm" 1 . Entangled ligh t exc ites a pure double-exciton 
state, J2fTf g (t)\f) [see eq. dA8k]. Consequently the purity 
of the density matrices in figure pK-B is tr/{^} = 1. The 
density matrix induced by stochastic light (figure [6p-D) is not 
in a pure state, we obtain tr/{^>j} « 0.24 independent of the 
pump frequency. Furthermore, it is apparent that by tuning 
the pump frequency, we can select the excitation of certain 
states. Figure [6}\ shows strong excitation of the states /n, 
/15 and fiQ, while in figure [6fe the states between fci and 




FIG. 8. (Color online) Second-order contribution to the total single- 
exciton manifold population g ee (t) , eq. |7J, with T — 30 fs. 

The total population of the single-exciton manifold J2 e Qee 
calculated to second order in the matter-field interaction [eq. 



(A3)] is displayed in figure [8] vs. the frequencies uj p and UJ2- 
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FIG. 9. (Color online) Total population of the double-exciton man- 
ifold J2f Qff [ ec l- Pi)] plotted vs. CJ2 and cj p . A: excitation by 
1 = 30 fs. B: Excitation by stochastic light. 



entangled light with . 



Due to the broadband nature of the beams, none of the peaks 
of the absorption spectrum in figure can be resolved. This 
population only depends on the spectral density, is identical 
for stochastic and entangled light, and is not susceptible for 
the non-classical features of entangled light. 
The total double-exciton population J2f Qff created by eq. 
( A8 ) is plotted in figure [9}\. The distribution of the single- 
exciton states to fourth order (see diagram II in figure|2]) looks 
almost identical, since eq. (TTOb includes an excitation to the 



double-exciton manifold, and is not shown. It depends only 
weakly on uj 2 , reflecting the energy uncertainty due to the time 
entanglement. Thus, the central frequency of each beam is not 
important, and in the following we set uo 2 = 11000 cm -1 . 
Two broad resonance bands, at 22000 — 23000 cm -1 (band 
a), and at 23000 - 24000 cm" 1 (band 6), dominate the dis- 
tribution. From our previous discussion of the density matrix, 
we anticipate that states /n.,/15 and f\§ dominate band a, 
whereas the states f 2 % — /27 seem to make up band b (see 
figure [6]). The signal induced by stochastic light (figure |9j3) 
shows no such structure. Due to the broad bandwidth, the 
level structure of the RC cannot be resolved at all. 

To get a more detailed picture of this state distribution, we 
depict the contribution of various double-exciton state popu- 
lations to the distribution [eq. (|A8|)] in figure [ToK. The dis- 



tributions are normalized at each uj p . States /n, /15 and fi$ 
dominate band a, whereas band b is dominated by states f 2 s, 
/ 2 4 and / 25 . At higher pump frequency, states f 27 , /30 and f 32 
are most pronounced. In general, the regions of leading contri- 
bution group around the states' energies, but the distributions 
may be asymmetric because of the presence of other states. 
We next turn to the distribution of single-exciton states (fig- 
ure [10^). Since the single excitons are obtained by emission 
from a two-exciton state, the distribution closely resembles 
the double-exciton distribution in figure [T0}\. For instance, 
when u p ~ 22200 cm -1 , the double-exciton state fu is on 
resonance. It decays primarily into single-exciton state e$, and 
accordingly the contribution of this state increases around the 
same value of the pump frequency. The same holds for, e.g., 
the double-exciton state /30, and the single-exciton state eio- 
However, in most cases the different contributions overlap, 
and the relation between the single-exciton and the double- 
exciton states is more complex. 
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FIG. 10. (Color online) A: contributions of the various f-states to the 
total population for CJ2 = 11000 era -1 with T — 30 fs upon 
excitation by entangled light. Insets: Excitation by stochastic light. 
B: same for the single-exciton states 



V. TWO-PHOTON-INDUCED FLUORESCENCE 





FIG. 11. (Color online) A: 2D spectrum of the fluorescence signal 
of entangled photons with entanglement time T — 30 fs. B: The 
corresponding spectrum with classical stochastic light. We set the 
frequency 002 — 11 000 cm -1 . Two-exciton peaks are indicated by 
vertical black lines, and single-exciton peaks by horizontal lines. 

Having discussed our simulations of the excited state dis- 
tribution created by entangled light, we now turn to the ques- 
tion whether it may be possible to detect fingerprints of these 
distributions in experiments. To that end, we calculate the 
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FIG. 12. (Color online) A: The dispersed fluorescence signal in- 
duced by entangled photons with entanglement time T — 30 fs 
and Up — 22160 cm -1 , B: enlarged region by the black box. C: 
Enlarged region for uj p — 22500 cm -1 . D & E: same as A & B, 
but for stochastic light. 




FIG. 13. (Color online) The fluorescence action spectrum (integrated 
over the emission frequency uo s ) vs. the pump frequency uj v \ signal 
created by entangled photons (solid, blue line), the stochastic light as 
(dashed, gray line), and signal created by two monochromatic beams 
at uj\ and CJ2 (dotted, red line). 

The double-exciton energies are plotted as solid, blue lines, and their 
heights are given by | J2 e l^gel^efl 2 • The single-exciton energies are 
indicated as dashed, red lines, with the height proportional to \fi ge | 2 . 



dispersed fluorescence signal. This involves transitions be- 
tween double- and single-exciton states as well as transitions 
between single-exciton states and the ground state, i.e. 

S f (uj s ;T) = ^2 \^fe\ 2 Pf(t',F)5(ujf -uj e - u 3 ), (23) 



S e (u s ; F) = ^ \/n eg \ 2 p e (t; T)S(u e - u s ). 



Here, u s is the emitted fluorescence frequency, and pf (p e ) the 
population of state f (e) given by eqs. ( [TT] ) and ( [TO} , respec- 
tively. More elaborate detection (gated time and frequency, 
photon statistics) is possible [49], but will not be considered 
here. 



In figure 11 we display this signal vs. the pump frequency 
uj p and the emission frequency uj s . The delta functions in 
( [23] ) and ( [24] ) were replaced by a Lorentzian, which means 
the width of the peaks in horizontal direction of figure [TT 
is instrumental. In analogy to the populations shown in fig- 
ure [9j the fluorescence simulation with entangled photons 
(fTT]A) shows two distinct resonances along the cj p -axis per- 
taining to bands a and b, respectively (see also figure 13). 
Along uj s (horizontal axis), it contains two contributions, one 
around uj s ~ 10000 cm -1 , and a higher-energy part between 
12000 and 13000 cm -1 , reminiscent of the absorption spec- 
trum. This pronounced structure can be exploited to enhance 
or suppress certain features, as will be demonstrated in the 
following. The simulation of the stochastic signal in figure 
[TTft shows no structure along the cj p -axis, and consequently 
does not allow a manipulation of the fluorescence signal. The 
dispersed fluorescence simulation for uj p = 22160 cm -1 
is shown in figure [T2]A. Most of the peaks can be directly 
attributed to the energy of single-exciton states, the transi- 
tion thus correspond to either fluorescence from the single- 
exciton manifold, or from double-exciton states, whose en- 



ergies closely match the sum of two single-exciton states. 
However, the simulations also show side peaks of the reso- 
nance at 10245 cm -1 , which cannot be assigned to single- 
exciton energies. The strongest of these peaks is highlighted 
by the black arrow in the inset figure [12^. It corresponds to 
the fn —> ej transition. Since state fu is most strongly ex- 



cited at the given pump frequency (see figure 10A), this peak 



is most pronounced. At pump frequency 22500 cm 



res- 



onantly excited states /15 and fiQ show different transitions 
into the single-exciton state e$ (see figure 
tion with stochastic light shown in figures 



12 



!). The simula- 
T2t) and E cannot 



resolve any of these sidepeaks due to its broadband nature. 



VI. CONCLUSIONS 

We had derived expressions for the single- and double- 
exciton density matrices of quantum systems interacting with 
arbitrary light sources perturbatively in the field. These were 
used to simulate the excitation of matter via entangled twin 
photons, and stochastic light with the same power spectrum. 
Applications to the reaction center of purple bacteria show 
that the populations strongly depend on the nature of the light, 
and we indicated how these properties could be observed in 
the frequency-resolved fluorescence measurements. The non- 
classical spectral profile of entangled light allows us to tar- 
get specific double-exciton states, and to explore all the ex- 
citation pathways to this state in a single shot. This shows 
up perspiciously in the fluorescence action spectrum in figure 



13 The spectrum reveals the level structure of the double- 
exciton manifold, whereas the stochastic light cannot resolve 
this structure. For completeness, we also plot the action spec- 
trum created by two cw laser beams with frequencies uj\ and 
co>2 at the peaks of the spectral density. Even though the narrow 
bandwidth of these beams allows for a good uj p -resolution, 
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this signal is most pronounced, when one of the laser beams 
is resonant with a single-exciton state, and thus closely resem- 
bles the linear spectrum in figure [5] 

The role of entanglement and coherent transport in the effi- 
ciency of exciton transport in photo synthetic complexes has 
been subject of intense debate I5QH571 . Using our formal- 
ism, one could also discuss the entanglement of quasiparticles 
in the double-exciton manifold (58l|59]|. This will require to 
decompose double-exciton states into single excitons, which 
goes beyond the scope of the present Paper. 
Other possible extensions are to excitation by pulsed light, and 
to include time- and frequency gating in order to obtain con- 
trol over the time and frequency resolution of the detection. 
The time resolution could allow us to monitor the dynamics of 
exciton transport and charge separation after interaction with 
entangled light in an experimentally feasible way. 
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Appendix A: Excitations induced by twin photons 

Using eq. ( [T2] ), the field correlation function of the leading- 
order contribution, given by eq. ([7]), yields 



(eHt 2 )E( Ti )) 



=Ctri 



T2 - Ti 



(n— t 2 ) _|_ e ~ iaj 2 (Ti— r 2 ) 



), (Al) 



where we have defined the triangular function 

\x\, 



tri(x) 



r 



o 



for \x\ < 1 
else. 



(A2) 



Expanding the matter correlation function into a sum-over- 
states expression, we obtain 



( 



+ sine ((w 2 - 0J ej g - ijej)T/2) \ , 



sine ((wi - uj eig + H e JT/2) 



(A3) 



where we dropped the constant factor C. The matrix element 
of the dipole operator connecting the states g and e is denoted 
H ge , the energy difference is given by u eg , and the lifetime 
broadening is denoted by j e . We can assume that r y ei w j ej . 
Since the entanglement time T is short with relative to other 
timescales in the system (T ~ 10 -3 cm -1 ), this contribution 
to the single-exciton states depends weakly on the pump fre- 
quency UJ p . 

The field correlation function of eq.([9]) can be recast as a nor- 
mally ordered term plus a commutator term, which we ne- 
glect. Using eq. ( p3] ) and dropping the constant factor C, we 
obtain for the populations 



2e,ei(*;r) =5R 



7e^ 4 



UJl - UJ e r g + lj e 



^2 ~ U eg + ij e 



-}■ 



(A4) 



This contribution is proportional to the entanglement time T, 
reflecting the fact that in ^ the time pairs (ti, ts) as well as 
(T2, T4) are correlated. This means that the entire process has 
to happen within the entanglement time T. Since this time is 
very short with respect to other timescales in the system, £ e?e j 
can be neglected. The only fourth-order contribution to the 
single-exciton manifold is thus given by ( [TO} . In a sum-over- 
states expression, it reads 



e',f a,b,c,a 



h 4 



^fej^gei 



^fei^gej 



J(uj a -u) eig +i~/ ei )T _ ^ e i(u b -u; e , g +i~f e ,)T 



Ul+^-UJfg+i^ff 0J ( 

-i(uj c -uj e g-ije )T _ 



I e -i(u d -uj e , g -ij e ,)T _ -|_" 



' lie 



(A5) 



where the summations a, b, c and d runs over uj\ and 002- It is peak is not. Additionally, taking into account 
apparent that the single-exciton peaks are multiplied by phase 
factors of the kind e *Oi-^e g +*7e)T _ w j 1 ji e ^ e two-exciton 

UJ ~ 0J eg + ij e 

2 iTe i^-^ g )T/2 sinc ^ _ Ueg )T/2) , (A6) 
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the single-exciton peaks have a very similar structure as ( |A3| ), 
and depend weakly on uj p . Thus, Q ei , ej (t) is dominated by the 
Lorentzian l/(u)i + 002 ~ ^fg + *7/)> an d this in turn means 
that spectroscopy with entangled photons can be used to probe 



the two-exciton manifold of aggregates. 
We further note that eq. (J - 
transition amplitudes Tf g ( 



We further note that eq. ( [TT] ) can be recast as the product of 



T fg (t) = j Jr fjr' (/(t)|0(r)0(rO|^)(O|^(r)^(rO|^) (A7) 



Clearly, the single-exciton resonances are suppressed, and 
the detuning between uo p and the two-exciton manifold com- 
pletely determines Qf.j (t;T). 

Appendix B: Electronic excitations with stochastic light 

Due to the different intensity scaling of classical signals, 
the main contribution to the signal will arise from the leading- 



order contribution in eq. ([7]). Since the stochastic light re- 
produces the spectral density of entangled light, we obtain the 
same result as for entangled light [see eq. ([A3}] . Additionally, 
one can filter out this contribution, and calculate ( [TQ| and ( [TT] ) 
to study the electronic excitations induced by higher-order in- 
teraction with stochastic light. 

Using eqs. ([9]), ( [TT] ) and ([22]), the fourth-order contributions 
to the single-exciton manifold yield 



£ ei , ej .n(t; r) = ^— ^2 /V'/V/sinc 2 ((u a - u e > g )T/2) sine 2 ((u b - u fe >)T/2) 



a,b e 



^gejVfa y^ Ufg _ Ue , g _ Ue _ g + ijf)(u fg - uj e 'g - uj ej9 - ijf) (uj e >e 3 - 2ij e )(u ejei - 2ij e ) 

+ /ifi ' ei/i/ej ((CJ /P - UJ e >g - UJ e%g ~ ilf){Ufg ^ g - UJ Cj g + ijf) (uJ e >e x ~ ! 

QfiJo r ) = VgeVef^ge>Ve>f 3 sinc 2 (K u eg )T/2) sine 2 ((uj b - (J fje )T/2) 



a,b e,e' 

1 1 

(U fjg ~ UJ eg ~ U e < Q - ilf){Uf ig ~ UJ eg ~ UJ e 'g + Hf) (^e'e + 2ij e )(u fifj + 2i<J f )) ' 



(B2) 



where we neglected the small imaginary part of the sinc 2 - 
functions. The summations a, b run over oj\ and U02, and we 
can assume that j e t w 7 e . In comparison with the excitations 
induced by entangled photons, we observe that the single- 
exciton peaks in (|A5|) and (|A8|) have a very similar structure 



to ( |B1| ) and ( |B2| ), respectively. However, eqs. ( [T0| ) and ( [TT] ) 
also contain Lorentzian factors that depend on uj\ +UJ2, which 
dominate the signal. Those resonances are replaced by a con- 
stant background contribution in the second lines of eqs. ( |B1| ) 
anddB2l). 
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